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Abstract 

Background: Argonaute (AGO) is the core component of RNA-induced silencing complex. The AGO gene family 
has been analyzed in various plant species; however, there is no report about AGOs in the well-known Traditional 
Chinese Medicine (TCM) plant, Salvia miltiorrhiza. 

Results: Through a genome-wide analysis, we identified ten SmAGO genes in 5. miltiorrhiza. Full-length cDNAs of all 
SmAGOs were subsequently cloned and sequenced. These SmAGOs were characterized using a comprehensive 
approach. Sequence features, gene structures and conserved domains were analyzed by the comparison of 
SmAGOs and AtAGOs. Phylogenetic relationships among AGO proteins from 5. miltiorrhiza, Arabidopsis and rice were 
revealed. The expression levels of SmAGO genes in various tissues of 5. miltiorrhiza were investigated. The results 
implied that some SmAGOs, such as SmAGOl, SmAG02, SmAG03, SmAG07 and SmAGOlO, probably played similar 
roles as their counterparts in Arabidopsis; whereas the others could be more species-specialized. It suggests the 
conservation and diversity of AGOs in plants. Additionally, we identified a total of 24 hairpin structures, representing 
six miRNA gene families, to be miRNA precursors. Using the modified 5 -RACE method, we confirmed that SmAGOl 
and SmAG02 were targeted by 5. miltiorrhiza miR168a/b and miR403, respectively. It suggests the conservation of 
AG01-miR168 and AGO2-miR403 regulatory modules in 5. miltiorrhiza and Arabidopsis. 

Conclusions: This is the first attempt to explore SmAGOs and miRNAs in 5. miltiorrhiza. The results provide useful 
information for further elucidation of gene silencing pathways in 5. miltiorrhiza. 



Background 

Small RNAs (sRNAs) involved in various gene silen- 
cing pathways play important and diverse roles in the 
development and differentiation of organisms through 
regulating gene expression at the transcriptional and 
post-transcriptional levels [1], affecting heterochromatin 
formation [2], and responding to biotic and abiotic stresses 
[3]. In plants, sRNAs are generated from double-stranded 
RNAs (dsRNAs) through various pathways and may be 
classified into two major classes, including microRNAs 



* Correspondence: sflu@implad.ac.cn 

Institute of Medicinal Plant Development, Chinese Academy of Medical 
Sciences & Peking Union Medical College, No.151, Malianwa North Road, 
Haidian District, Beijing 100193, China 



(miRNAs) and small interfering RNAs (siRNAs), based on 
the source of dsRNA [4]. miRNAs are produced from plant 
transcripts with internal stem-loop structures, whereas 
siRNAs are derived from transcripts with inverted-repeat 
sequence, dsRNAs copied from single-stranded RNA 
(ssRNA), over-lapping regions of bidirectional transcripts, or 
dsRNAs formed by virus replication. The biogenesis path- 
ways of plant sRNAs involve in various gene families, such 
as the Dicer-like (DCL) family and the RNA dependent 
RNA polymerase (RDR) family, and each pathway appears 
to be taken part in by different member of a gene family [4]. 
To regulate gene expression, the generated sRNA duplexes 
from dsRNAs are loaded into RNA-induced silencing 
complexes (RISCs) with Argonautes (AGOs) as the central 



O© 2013 Shao and Lu; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
BiolVlGCl C6ntTcll Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly cited. 



Shao and Lu BMC Genomics 2013, 14:512 
http://www.biomedcentral.eom/1 471 -21 64/1 4/512 



Page 2 of 12 



components [5]. RISCs remove the star strand (known as 
miRNA* or siRNA*) of sRNA duplex and select the func- 
tional strand as a guide to interact with homologous RNA 
or DNA molecules for direct RNA cleavage, translational 
repression or DNA methylation [6,7]. 

AGOs are ribonucleases with two conserved domains, 
including PAZ and PIWI [8]. The PAZ domain contains 
a specific binding pocket that can anchor sRNA du- 
plexes with two-nucleotide 3' overhang. The PIWI do- 
main exhibits endonuclease activity and the structure of 
PIWI domain folded is similar to RNase H [9]. The 
endonuclease activity of PIWI domain is performed by 
an active site usually carrying an Asp-Asp-His (DDH) 
or Asp-Asp-Asp (DDD) motif. AGOs are usually 
encoded by a multiple gene family in organisms and the 
number of AGO genes differs in different organisms 
[10]. Fission yeast has only one AGO gene, whereas 
insects, mammals and worms have five, eight, and 
twenty six AGO genes, respectively. Annotation of the 
Arabidopsis and rice genomes revealed ten and eighteen 
AGO genes, respectively [11]. Additionally, a total of 
eighteen maize and fifteen tomato AGO genes have been 
identified [12,13]. Although a large number of AGO 
genes have been found in various plants, the majority 
were predicted by computational approaches based on 
sequence similarity. Only a small proportion of known 
AGOs were identified or confirmed by full-length cDNA 
cloning. Among them, Arabidopsis AGOl is the most 
well-studied plant AGO gene. It encodes the core com- 
ponent of RISCs associated with the action of miRNAs, 
trans-acting siRNAs (ta-siRNAs) and transgene-derived 
siRNAs [14]. AG02 protein is involved in antiviral 
defense by catalyzing viral RNA cleavage in Arabidopsis 
plants [15]. Arabidopsis AG04, AG06 and AG09 genes 
recruit endogenous 24nt sRNAs for DNA methylation, 
which causes target gene silencing at the transcriptional 
level [16,17]. AG07 is involved in the generation of 
ta-siRNAs from TAS3 by collaborating with miR390 in 
Arabidopsis [18]. Arabidopsis AGO10 gene modulates 
shoot apical meristem maintenance and establishment of 
leaf polarity by repressing miR165/166 [19,20]. In 
addition, although AGO proteins are central compo- 
nents of RISCs involved in sRNA-mediated RNA cleav- 
age, translational repression and DNA methylation, 
some of them are regulated by sRNAs through the feed- 
back mechanism. For instance, the expression level of 
AtAGOl is regulated by miR168 through direct cleavage 
of AtAGOl transcripts [21]. Similarly, AtAG02 is regu- 
lated by miR403 [22]. The function and regulatory 
mechanism of AGO genes from plant species other than 
Arabidopsis is largely unknown. 

S. miltiorrhiza, which produces two major groups of 
bioactive compounds, lipophilic diterpenoid tanshinones 
and hydrophilic phenolic acids, is a well-known tradi- 



tional Chinese medicine (TCM) widely used for treating 
various human diseases, such as dysmenorrhoea, amen- 
orrhoea and cardiovascular disease, for thousands of 
years [23,24]. It is also an emerging model plant for 
TCM studies because of its relatively small genome size, 
short life cycle, undemanding growth requirements, and 
significant medicinal value [25,26]. The S. miltiorrhiza 
genome has been preliminarily decoded and a working 
draft of the genome is currently available (Chen et al, 
unpublished data). The interest of gene silencing path- 
ways in S. miltiorrhiza is increasing. With the aim to 
elucidate the core components of gene silencing path- 
ways, we performed a genome-wide prediction of the S. 
miltiorrhiza AGO gene family. Molecular cloning of 
AGO genes was carried out for validation and error cor- 
rection of computational prediction. The characteristics 
of S. miltiorrhiza AGOs were revealed by a comprehen- 
sive analysis, including comparison with AGOs from 
other plant species, gene expression profiling, and ana- 
lysis of posttranscriptional regulation. The results pro- 
vide useful information for further elucidation of gene 
silencing pathways in S. miltiorrhiza. 

Results 

Genome-wide prediction of S. miltiorrhiza AGO genes 

To predict S. miltiorrhiza AGO genes at the genome level, 
we downloaded all of 10 Arabidopsis and 19 rice AGO 
amino acid sequences from GenBank (http://www.ncbi. 
nlm.nih.gov/protein). BLAST analysis of Arabidopsis and 
rice AGOs against the current assembly of the S. 
miltiorrhiza genome (Chen et al, unpublished) was then 
performed using the tBLASTn algorithm [27]. An e-value 
cut-off of 10" 10 was applied to the homologue recognition. 
As a result, a total of 10 genomic loci of SmAGO genes 
were identified. The genomic DNA sequence was re- 
trieved and gene models of 10 SmAGO genes were pre- 
dicted using Genscan (http://genes.mit.edu/GENSCAN. 
html) [28]. The models were further examined and 
corrected manually by comparison with AGO genes iden- 
tified from other plant species using the BLASTx algo- 
rithm (http://www.ncbi.nlm.nih.gov/BLAST) [27]. All of 
10 deduced proteins share high sequence similarity with 
known plant AGOs and contain the conserved PAZ and 
PIWI domains, suggesting they are authentic AGOs. The 
identified AGO genes are named SmAGOl to SmAGOlO, 
respectively, based on high sequence similarity with corre- 
sponding Arabidopsis AGOs. 

Sequence feature, gene structure and conserved domain 
comparison of SmAGOs and AtAGOs 

It is very important to know the correct cDNA sequence 
for systematic characterization of AGOs in S. miltiorrhiza. 
In order to confirm the results from prediction and 
correct errors of computation, molecular cloning of full- 
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length SmAGO cDNA was carried out. As a result, all of 
ten predicted SmAGO genes were experimentally vali- 
dated. Analysis of sequence features showed that the 
length of open reading frames (ORFs) of SmAGOs varied 
from 2,472 (SmAGOS) to 3,195 bp (SmAGOl) (Table 1). 
The length of 5 ' and 3 ' UTRs was between 30 and 202 bp 
and between 91 and 337 bp, respectively. The size of de- 
duced SmAGO proteins varied between 823 and 1064 
amino acids, the molecular weight (Mw) varied from 93.2 
to 118.7 kDa, and the theoretical p/was between 9 and 10 
(Table 1). These sequence features are quite similar to that 
of AtAGOs in Arabidopsis (Table 2). The similarity be- 
tween SmAGOs and AtAGOs was also observed in gene 
structures. Many SmAGOs genes have the number and 
phase of introns similar to an AtAGO, such as SmAGOl/ 
AtAGOl, SmA G04/AtA G04, SmAG07/AtAG07, SmA 
G08/AtAG08, and SmAG09IAtAG09 (Tables 1 and 2, 
Figures 1 and 2). Additionally, both SmAG06/AtAG06 
and SmAG09IAtAG09 have an intron in 3' UTR 
(Figures 1 and 2). The results suggest the conservation 
between SmAGOs and AtAGOs, 

Search S. miltiorrhiza AGO proteins for conserved do- 
mains against the NCBI Conserved Domain Database 
(CCD) revealed that all SmAGOs contained DUF1785, 
PAZ, and PIWI domains (see Additional file 1), which 
were consistent with the results from other plant AGOs 
[29]. The function of DUF1785 is currently unknown; 
whereas PAZ has been shown to play roles in binding 
sRNA duplexes and PIWI are important in cleavage of 
target RNA [30,31]. The endonuclease activity of PIWI 
domain is performed by an active site usually carrying a 
DDH or DDD motif [15,32]. Additionally, a conserved 
histidine at position 798 of Arabidopsis AGOl was 
found to be critical for in vitro endonuclease activity 
[14]. Analysis of ten S. miltiorrhiza AGOs showed that 
half of them, namely SmAGOl, SmAG02, SmAG03, 
SmAG07 and SmAGO 10, contained the conserved 
DDH/H798 or DDD/H798 residues; whereas the con- 
served residues were not observed in the other SmAGOs 



(see Additional file 2). In SmAG04, SmAG05 and 
SmAG08, the third histidine was missing or replaced by 
leucine. SmAG06 and SmAG09 possessed the conserved 
DDH triad but histidine at 798th position was either re- 
placed by alanine or proline (see Additional file 2). 

Phylogenetic analysis of AGO proteins in S. miltiorrhiza, 
Arabidopsis and rice 

Phylogenetic analysis using the PAZ and PIWI domains 
for rice, Arabidopsis, Caenorhabditis elegans, Drosophila 
melanogaster, and mouse had previously revealed that 
animal AGOs clustered into two subgroups: Al and A2, 
whereas all plant AGOs could be divided into four sub- 
groups: AGOl, ZIPPY, AG04, and MEL1 [33]. The 
AGOl and MEL1 subgroups had a common lineage with 
Al, whereas plant ZIPPY and AG04 subgroups and ani- 
mal A2 subgroup could be diverged from an ancestral 
lineage [11]. In order to determine the evolutionary rela- 
tionship of S. miltiorrhiza AGOs, full-length AGO pro- 
tein sequences from S. miltiorrhiza, Arabidopsis and rice 
were aligned and an unrooted neighbor-joining tree was 
constructed. The results showed that ten SmAGOs 
could also be divided into four subgroups (Figure 3). 
Moreover, the clustering remained similar when only the 
PAZ and PIWI protein domains were used for phylogen- 
etic analysis (data not shown). SmAGOl and SmAGO 10 
are included in the AGOl subgroup with Arabidopsis 
AtAGOl and AtAGOlO, and rice OsAGOla-OsAGOld 
and OsPNHl. SmAGOl shares high similarity with 
AtAGOl associated with the action of miRNAs, ta- 
siRNAs and transgene-derived siRNAs [14], whereas 
SmAGO 10 is highly similar to AtAGOlO modulating 
shoot apical meristem maintenance and establishment of 
leaf polarity [20]. SmAG02, SmAG03 and SmAG07 are 
members of the subgroup ZIPPY. Similarly, three 
Arabidopsis AGOs, including antiviral defense-associated 
AtAG02, function-unknown AtAG03 and AtAG07 in- 
volved in the generation of ta-siRNAs from TAS3 [15], are 
also included in this subgroup. SmAG07 shares greater 



Table 1 Sequence features and intron number of SmAGOs 



Gene name 


Accession number 


cDNA (bp) 


ORF (bp) 


5'UTR (bp) 


3'UTR (bp) 


Protein (aa) 


Mw (kDa) 


p/ 


Intron no. 


SmAGOl 


KF1 53679 


3429 


3195 


121 


113 


1064 


118.7 


9.68 


20 


SmAG02 


KF1 53680 


3369 


2964 


68 


337 


988 


110.1 


9.85 


1 


SmAG03 


KF1 53681 


2975 


2793 


86 


96 


930 


106.0 


9.86 


1 


SmAG04 


KF1 53682 


2974 


2613 


81 


280 


870 


97.7 


9.16 


21 


SmAG05 


KF1 53683 


2905 


2472 


177 


256 


823 


93.2 


9.64 


20 


SmAG06 


KF1 53684 


3212 


2712 


172 


328 


903 


101.1 


9.37 


12 


SmAG07 


KF1 53685 


3182 


3045 


46 


91 


1014 


115.5 


9.28 


2 


SmAG08 


KF1 53686 


3075 


2673 


202 


200 


890 


100.0 


9.79 


18 


SmAG09 


KF1 53687 


2915 


2667 


30 


218 


888 


98.7 


9.24 


21 


SmAGOl 0 


KF1 53688 


3118 


2907 


39 


172 


968 


108.9 


9.44 


18 
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Table 2 Sequence features and intron number of AtAGOs 



Gene name 


Gene model 


cDNA (bp) 


ORF (bp) 


5'UTR (bp) 


3'UTR (bp) 


Protein (aa) 


Mw (kDa) 


p/ 


Intron no. 


AtAGOl 


AT1G48410.1 


3601 


3147 


267 


187 


1048 


116.2 


9.65 


21 


AXAG02 


AT1G3 1280.1 


3343 


3045 


66 


232 


1014 


113.4 


9.97 


2 


AtAG03 


AT1G3 1290.1 


3585 


3585 






1 194 


129.2 


9.78 


2 


AtAG04 


AT2G27040.1 


3170 


2775 


134 


261 


924 


102.9 


9.19 


22 


AtAG05 


AT2G27880.1 


2994 


2994 






997 


111.1 


9.92 


19 


AtAG06 


AT2G32940.1 


2997 


2637 


139 


221 


878 


98.7 


8.59 


22 


AXAG07 


AT1G69440.1 


3160 


2973 


150 


37 


990 


113.4 


9.59 


2 


AtAG08 


AT5G21 030.1 


2553 


2553 






850 


95.5 


8.92 


20 


AXAG09 


AT5G21 150.1 


3017 


2691 


115 


211 


896 


100.5 


9.45 


22 


AtAGOl 0 


AT5G43810.1 


3330 


2967 


128 


235 


988 


110.9 


9.64 


18 



similarity with AtAG07. SmAG02 and SmAG03 are 
highly similar to AtAG02 and AtAG03. The MEL1 sub- 
group contains only one S. miltiorrhiza AGO, SmAG05, 
which is similar to the function-unknown Arabidopsis 
AtAG05. On the contrary, the AG04 subgroup is the big- 
gest among four plant AGO subgroups. It contains five S. 
miltiorrhiza AGOs, including SmAG04, SmAG06, 
SmAG08 and SmAG09. 

Differential expression of SmAGO genes 

The expression of SmAGO genes in flowers, leaves, 
stems and roots of 2-year-old, field nursery-grown S. 
miltiorrhiza was analyzed using quantitative RT-PCR 
technology. SmllBQlO was chosen as an endogenous 
control as previously described [26]. All of ten SmAGOs 
were expressed in S. miltiorrhiza tissues analyzed, 
whereas differential expression patterns were observed 
(Figure 4). SmAGO 1 and SmAGO 10 clustered in the 
AGOl subgroup were expressed in flowers, leaves, stems 
and roots as their Arabidopsis counterparts, AtAGOl 
and At AGO 10, respectively [11]. It is consistent with the 
ubiquitous roles of AGOs in the AGOl subgroup. 
Although SmAG02, SmAG03 and SmAG07 were clus- 
tered in the ZIPPY subgroup, their expression patterns 



were distinct (Figure 4). On the other hand, the expres- 
sion patterns of SmAG02, SmAG03 and SmAG07 were 
quite similar to their Arabidopsis counterparts, AtAG03, 
AtAG02, and AtAG07, respectively [11]. By contrast, 
the expression pattern of SmAGO S was distinct with 
their Arabidopsis and rice counterparts in the MEL1 
subgroup. SmAGOS showed the highest expression in 
roots, followed by stems, less in flowers and leaves. It is 
distinct with its counterparts, AtAGOS in Arabidopsis 
and OsMELl, OsAG012, OsAG013 and OsAG014 in 
rice, showing specifically expression in reproductive tis- 
sues (Figure 4) [34]. Similar to the SmAGOs in the 
ZIPPY subgroup, SmAG04, SmAG06, SmAGOS and 
SmAG09 clustered in the AG04 subgroup exhibited di- 
vergent expression patterns (Figure 4). SmAGOS was 
more flower-specific, an expression pattern similar with 
SmAG07 in the ZIPPY subgroup; whereas, SmAG06 
and SmAG09 were more root-specific. The results indi- 
cate the functional conservation and diversity of AGOs. 

miRNA-mediated posttranscriptional regulation of SmAGO 
genes 

To determine whether S. miltiorrhiza AGOs are regu- 
lated by miRNAs, we performed a target search of plant 
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Figure 1 Gene structures of AGOs in S. miltiorrhiza. Exons, UTRs, introns and intron phases are shown. 
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exon UTR intron 0, 1, 2: intron phase 

Figure 2 Gene structures of AGOs in Arabidopsis. Exons, UTRs, introns and intron phases are shown. 



miRNAs in miRBase against ten full-length SmAGO 
cDNA sequences using psRNATarget [35,36]. The 
maximum expectation of 3.5 was applied in the target 
search. A total of 31 miRNA families were found to have 
perfect or near-perfect complementarity to SmAGOs. It 
includes 28 mature miRNAs and 3 likely miRNA* (aly- 
miR167d-3p, ptc-miR169n-3p, and gma-miR396j). Plant 
miRNA sequences belonging to the 31 families were 
aligned with the current assembly of the S. miltiorrhiza 
genome using SOAP2 with two mismatches allowed 



[37]. Genomic DNA fragments surrounding these 
miRNA sequences were used to predict the secondary 
structure using the mfold program [38]. A total of 24 
hairpin structures were identified for six miRNA 
families, including miR167, miR168, miR169, miR396, 
miR403 and miR530 (Figure 5). Manual examination of 
the complementarities between SmAGOs and the identi- 
fied S. miltiorrhiza miRNAs/miRNAs* and calculation of 
penalty scores as previously described [39] showed that 
the scores for S. miltiorrhiza miR168:«SmAGOi, miR403: 
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Figure 4 Expression of SmAGOs in flowers (Fl), leaves (Le), stems (St) and roots (Rt) of S. miltiorrhiza. Fold changes of SmAGO expression 
are shown. Expression levels were quantified by qRT-PCR. SmUBQIO was used as a reference gene. The levels in roots were arbitrarily set to 1. 
Error bars represent the standard deviations of three technical PCR replicates. 



SmAG02 and miR530:SmA GOl were 2, 0, and 3.5, re- 
spectively. However, the scores for miR167*:«SmA G08, 
miR169*:SraAG05, miR396*:SraA G03 and were at least 
5.5, 9.5, and 5.0, respectively. High penalty scores for S. 
miltiorrhiza miRN A* \SmAGOs were due to low conser- 
vation among plant miRNAs*. To verify whether these 
miRNAs can mediate the cleavage of SmAGO tran- 
scripts, we isolated RNAs from roots and pooled sam- 
ples containing flowers, leaves, stem and roots of S. 
miltiorrhiza and performed the modified 5 '-rapid ampli- 
fication of cDNA ends (RACE) for SmAGOs. The 5'- 
RACE products revealed that SmAGOl and SmAG02 
are indeed the targets of S. miltiorrhiza miR168 and 
miR403, respectively (Figure 6). miR168 targets to the 
DUF1785 domain (from 477 to 497 bp) region of 
SmAGOl, while the target site of miR403 is within 
the 3'UTR of SmAG02. It is consistent with the re- 
sults from Arabidopsis [21,22], suggesting the conser- 
vation of miR168- and miR403-mediated regulation of 
AGOs between S. miltiorrhiza and Arabidopsis. We 
also cloned 5 '-RACE products with the 5 '-end 
mapped to the miR169* complementary region of 
SmAGOS; however, the positions are not between the 
10 and 11 nucleotides from the 5 '-end of the miRNA 
(data not shown). No 5 '-RACE products were 
obtained for SmAG08, SmAGOS and SmAGOl that 
were predicted to be targeted by miR167, miR396, 



and miR530, respectively. In addition, miRNA-specific 
qRT-PCR [40] was performed to examine the expres- 
sion patterns of S. miltiorrhiza miR168 and miR403 
experimentally validated to target SmAGOl and 
SmAG02 for cleavage. The results showed that 
miR168 and miR403 were expressed in all tissues an- 
alyzed. The levels of miR168 in roots and flowers 
were about three times the levels in leaves and stems. 
The level of miR403 in flowers was more than twice 
the level in other tissues analyzed (Figure 7). 

Discussion 

Identification of ten full-length AGO cDNAs in 
S. miltiorrhiza 

Although Argonautes play very important roles in 
small RNA-mediated gene silencing and a large num- 
ber of AGO genes have been found in plants, many 
of them were identified through computational pre- 
diction based on sequence similarity. For example, 
often Arabidopsis AtAGOs, seven, including AtAGOl, 
AtAG02, AtAG04, AtAG06, AtAG07, AtAG09 and 
AtAGOl 0, have been experimentally tested, whereas 
the other three, including AgAG03, At AGOS and 
AtAG08, were predicted computationally (http://www. 
arabidopsis.org/). Among nineteen rice OsAGO genes, 
only six, including OsAGO la, OsAGO lb, OsAGO lc, 
OsAGOld, OsAG07 and OsPNHl, have been cloned 
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MIR167a (AG=-46.0) 

CUU-- U UAU uuccuuu 

5' -GAGAGGCUGAAGCUGCCAGCAUGAUCUG UCCA GUUUUUUAAA UUAAUCCA U 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I II I I I I I I I I 

3'-CUCACUAACUUCGACGAUCUUACUGGAC GGGU UAAAAAGCUU AAUUAGGU / 

UACGC - UU- UCUUCAC 



MIR167d (AG=-43.5) 

G UAA AGG 

5' -UAGUAGUUGAAGCUGCCA CAUGAUCU CUUGCCUC C 



MIR167b(AG= -52.1) 

AUA GUUUCUUCUUCGCUUUCUUCUAGAAAGA 
5' -GAAUAAGUGAAGCUGCC AGCAUGAUCUACGUCGGG GUUU A 



MIR167e (AG= -52.7) 

UGUC G AG 

5'- GGGUAAAGCUGCCA CAUGAUCUAAGCUUUCCUUCCU U 



MIR167C (AG= -52.1) 

CU CAACAAC UCC UUAAACGCUUCAUAAAUCUGGCCAUUUUGAAUCCCAUUUUCUU 

5' -GAGAGGUU AAGCUGCCAG UGAUCUGAUUA ACC AUUUUUGAGAOU C 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I III I I I I I I I I I I I I 

3'-CUCACUAA UUCGACGGUC ACUGGACUAGU 

C U- CACCU- 



MIR168a (AG= -61.4) 

U GUUGAAUCUCUACGUAUGCGUGUGAGUAUA 
5' -CUCUAGUUCGCUUGGUGCAGGUCGGGAACUU AUAUUC GC U 



MIR169a (AG= -47.3) 

U CG- 
5'- GGAUUGCAGCCAAGGAUGACUUGCCGGC 



-GAGGCUAAGUCAACUACGUUCCGCCCUUGGA UAUAAG CG 



GUUUAUGUCGCCGUUAGUUUACAAUGCGGC 



MIR168b (AG= -60.6) 

5 ' -CUCUAAUUCGCUUGGUGCAGGUCGGGAACUU 



UCUUCGCCGCCGCCGCGCUCCAAUCACCAAUUGCGUGUUU 



3 ' -GAGGCUAAGUCAACUACGUUCCGCCCUUGAA 



CUAAAAAGGUGGCUCAGCCUUAGCUCCUAACGGCGUCGUG 



MIR169e (AG= -39 5) 



AUGUAAUUAUAU 



AGAUGUAGCCAAGGAUGACUUGCCGG 



UCCGUAUCGGUUUUUAUUGAACGGCC 



AAGGUUUGUUAC 



MIR169b (AG= -49.6) 

U U UG- ACUC CA 

5'-AGUUUGA UAGCCAAGGA GACU GCCUG CCAA AUCAAAACUCGAAAC U 

3'-UCACUCU AUCGGUUCCU CUGA CGGAU GGUU UAGUUUUUGGUUUUG C 
A - - UUA GUUU UA 



MIR169C (AG= -52.2) 

U - U U CU- UUG UA 

5'- GUUUGG UAGCCAAGGA GACU GCCUG CCAUUGA AUUAAUUCCAUAUAUAUA U 

3'- CAGUUC AUCGAUUCCU CUGA CGGGU GGUGGUU UAAUUAGGGUAUAUAUGU G 
A A - - UUU — UU 



MIR169d (AG= -45.7) MIR169g (AG= -36.2) 

GGUUU - U U ACCCUCCAUUAUUGGGGGUUUUCUAGACUUUUACU C UAG- A U CC CA 

5'- GA UAGCCAAGGA GAC UGCCUGG U 5'- GUU UAGCCAAGA UGACU GCCCCAAU AUAUGUG 

I I I I I I I I I I I I I M I I I I I I I III I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3' - CU AUCGGUUCCU CUG GCGGACC / 

CACU A - - ACAGGUUUAACCUAUAAAACCCACAGUUUUAUAAC 



MIR169f (AG=-40.9) 

U G - UUAUUGAGAUAAUAAGAAGGCCUCACAAUAAAUUUGUGCUAUCUAUCUUUCUAUCUAUCUAUCUAUCU 

5'- GUUUAUU AGCCAAGAAUGACUU CCUC AAAUC A 



■ CAGUUGA UCGGUUCUAGCUGAA GGAG UUUAG 



U 



MIR169h (AG= -35.1) 

G U U U UAAU 

5'- GUU AGUAGCCAAGAA GAUU GCCUGGAAAUCUA I 



UUCUCGCACUCCGUACGUUGUUUAUUUAUCUAUUUUACUAUCAUAUUAUAUACCUGUGUCUAUCUAUC 
MIR169I (AG= -34.1) 

AG- UU AU 

5'-AAUUUGG UAGCCA GAUGAC GCCUAAAUAUUUUACAGGUUUCAUA 



III 



: I 



I 



I 



3' -UUAAACU AUCGGU CUACUG CGGAUUUUUGGUAUAAUCAAAAUAU 
U CUA -- Al 



MIR169j(AG=-43.6) 

C - - UU -- A- UAU CA 

5'- AUUUG GUAGCCAA GGAUGAC GCCUA AAU UGAUU GGCUUCUCAU A 



MIR169k (AG= -31.3) 

AU - CUAU-- 

5'- UUGGUUAGCCAAGAAUGA CUUGC 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3'- AGUUUAUCGGUUCCUGCU GAACG UGGUUUGUUUUG 
U C AAUCACUU Ui 



MIR396a (AG= -47.8) 

U ACAAAUU- C 

5'- CAUGCUUUUCCACAGCUUUCUUGAACU GCUUCGCUUCUUCAU G 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I II 

3'- GUACAAAAGGGUGUCGAAAGAACUUGA UGGGGUUAAGAACUA A 

G AUAUAUUU G 



MIR396b (AG= -50.0) 

CU C- A - C 

5'- UAUUCUUCCACAGCUUUCUUGAACUGCG CGA UUGUA CGUCGCCG A 

III I I II I I I I I I I I I I I I I I I I I I III I I I I I I I I I III 
3'- AUAGAAGGGUGUCGAAAUAACUUGGCGU GCU AACGU GCGGAGGC / 

GC AU - A C 



MIR396C (AG= -39.8) 

CU - UGCC 

5' - UAUUCUUCCACAGCUUUCUUGAACUGCA GAUC ACUGUA U 

3'- AUAGAAGGGUGUUGAAAUAACUUGAUGU CUAG UGACGU G 

GC UUCUUA U CAUC 



MIR396d (AG= -42.2) 

UUUUAUGAUGAAUAUUAU 
5' -CAUGCUUUUCCACAGCUUUCUUGAACUUC U 

3' -GUACAGAAGGGUGUCGAAAGAACUUGAAG / 
CAGUUCCAUUUUAAGAAC 



MIR403 (AG= -49.4) 

UCAU UCAU-- 
5'-CAUAUU GUUUGUUCGUGAGUCUAACGU 

3'-GUAUAG CAAA 
UGCU 

MIR530 (AG= -73.8) 

GC — AAUUUAGUAUUUGUUUAAACUUUUAUUAUGAAUUUUUUCUUUGUUUUUUCGUCAAUGCUUCGUUUCAUUCACCAGCA 
S'-CUUUACCUGCAUUUGCACCUGCACCUUCUU UUCUUCGC C 

3' -GAAGUGGACGUAUACGUAGACGUGGAGGAA AAGAGGCG / 
GUUA AAGGUACAAUUGCUCGAGUCAUUAACACACACACACACACACAAGAUCUGAAGAUCGUCGUAAUGUGAAUUUUAAGCA 

Figure 5 Predicted hairpin structures of S. miltiorrhiza miRNA precursors. Mature miRNA sequences are indicated in red. miRNA* are 
indicated in blue. The red, green and blue vertical lines indicate G:C, G:U and A:U pairing, respectively. 



(http://www.ricedata.cn/gene). Except for SlAGOl-1, 
SlAGOl-2 and SIAG07, twelve of fifteen tomato 
SIAGO genes have not been experimentally confirmed 
[41,42]. In this study, we performed a genome-wide 
prediction of ten SmAGOs using computational ap- 
proaches, and then cloned the full-length cDNAs of 
all predicted SmAGOs, The number of identified S. 
miltiorrhiza AGO genes is comparable with that in 
Arabidopsis, although it is significantly less than the 



number in rice and maize, which are nineteen and 
eighteen, respectively [11,12]. It indicates that less 
duplication events are occurred for AGO genes in S. 
miltiorrhiza and Arabidopsis as compared with rice 
and maize AGOs, most of which are evolved by 
duplication events [11,12]. The results provide very 
useful information for further elucidation of AGO 
functions in S. miltiorrhiza and gene model predic- 
tion of AGOs in other plant species. 
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...GUCCCGAGCUGCACCAAGCUA... ...CGAGUUUGUGCGUGAAUCUAA... 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

3' -AAGGGCUGGACGUGGUUCGCU-5' miR168a/b 3 ' -GCUCAAACACGCACUUAGAUU-5 ' miR403 

Figure 6 Experimental validation of miR168a/b- and miR403-mediated cleavage of Sm AGO! and SmAG02, respectively. Cleavage sites 
were determined by the modified 5' RNA ligase-mediated RACE. Heavy black lines represent ORFs. The lines flanking gray regions represent 
nontranslated regions. miRNA complementary sites with the nucleotide positions of Sm AG0 1 and SmAG02 cDNAs are indicated. The RNA 
sequence of each complementary site from 5' to 3' and the predicted miRNA sequence from 3' to 5' are shown in the expanded regions. 
Watson-Crick pairing is indicated by vertical dashes. Vertical arrows indicate the 5' termini of miRNA-guided cleavage products, as identified by 
5 -RACE, with the frequency of clones shown. 



Conservation and diversity of SmAGOs and AtAGOs 

Plant AGO proteins share three highly conserved do- 
mains, including DUF1785, PAZ, and PIWI [43]. Con- 
sistently, all SmAGOs were found to contain these 
domains (see Additional file 1). PAZ functions in bind- 
ing sRNA duplexes and PIWI is involved in RNA cleav- 
age, whereas the function of DUF1785 remains to be 
elucidated [30,31]. The conserved DDH/H798 or DDD/ 
H798 residues in PIWI domain have been demonstrated 
to be critical for the endonuclease activity of AGO pro- 
teins [15,32]. The conserved residues were also found in 
five of ten SmAGOs (see Additional file 2). It includes 
SmAGOl and SmAGOlO belonging to the AGOl sub- 
group and SmAG02, SmAG03 and SmAG07 included 
in the ZIPPY subgroup. Consistently, expression profil- 
ing of these SmAGOs revealed similar patterns with their 
Arabidopsis counterparts (Figure 4) [11]. It indicates S. 
miltiorrhiza AGOs in the AGOl and ZIPPY subgroup 
may play similar functions in the action of miRNAs, 



ta-siRNAs and transgene-derived siRNAs and in antiviral 
defense as their counterparts in Arabidopsis [15]. 

On the other hand, in the other five SmAGOs, one or 
two of the conserved residues was missing or replaced 
by other residues. It includes SmAG05 belonging to the 
MEL1 subgroup, and SmAG04, SmAG06, SmAG08 
and SmAG09 included in the AG04 subgroup. Com- 
parison of AGOs revealed that many Arabidopsis and 
rice AGOs belonging to the MEL1 and AG04 subgroups 
were also lack of the conserved DDH/H798 or DDD/ 
H798 motif, such as Arabidopsis AtAG04, AtAG06, 
AtAG08 and AtAG09, rice OsAG04a, OsAG04b, 
OsAG015 and OsAG016 included in the AG04 sub- 
group, and rice OsAGOll, OsAG013, OsAG014 and 
OsAG018 belonging to the MEL1 subgroup [11,13]. 
The function of AGOs in the MEL1 subgroup is cur- 
rently unknown, whereas three of four Arabidopsis 
AGOs included in the AG04 subgroups have been 
shown to recruit endogenous 24nt sRNAs for DNA 
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methylation [44-46]. It indicates that some of SmAGOs 
in the MEL1 and AG04 subgroups probably play a role 
in DNA methylation instead of RNA cleavage in S. 
miltiorrhiza. Since the expression patterns of SmAGOs 
in the MEL1 and AG04 subgroups are distinct with 
their Arabidopsis counterparts (Figure 4) [11], some 
SmAGOs in these subgroups may play more species- 
specialized roles. Further elucidation of these species- 
specialized roles will definitely add new insights into 
AGO-associated gene silencing. In addition, absence of 
the conserved DDH/H798 or DDD/H798 motif was also 
found in some AGO proteins from other plant species, 
such as maize and tomato [12,13]. It is possible that 
some of the AGO proteins without the deeply conserved 
DDH/H798 or DDD/H798 motif are still capable of tar- 
get RNA cleavage. In this case, the unconserved residues 
in the motif may contribute to endonuclease activity. 
Otherwise, they may be not a cleavage component if the 
deletion or replacement of conserved residues in the 
motif results in loss of endonuclease activity. Further ex- 
perimental analysis is necessary to clarify the hypothesis. 

S. miltiorrhiza miRNA identification and 
posttranscriptional regulation of SmAGO genes 

miRNAs are a class of small endogenous non-coding 
RNAs with size about 21 nucleotides. They are derived 
from primary miRNAs (pri-miRNAs) transcribed from 
miRNA loci [47]. Pri-miRNAs have internal stem-loop 
structures that are cleaved by DCL1 to form miRNA 
precursors, known as pre-miRNAs. miRNAs play vital 
roles in plant development and stress responses 
and have been identified from various plant species 
[39,48,49]; however, there is no report for miRNAs in S. 
miltiorrhiza, an emerging model medicinal plant. 
Through a computational approach, we predicted 24 
miRNA genes in S. miltiorrhiza for the first time. They 
represent 6 miRNA gene families, including miR167, 
miR168, miR169, miR396, miR403 and miR530 (Figure 5). 
Using the modified 5 '-RACE method, we confirmed 
that SmAGO 1 and SmAG02 were targeted by S. 
miltiorrhiza miR168a/b and miR403, respectively 
(Figure 6). Consistently, in Arabidopsis, AtAGOl and 
AtAG02 were also regulated by miR168 and miR403, 
respectively [21,22]. It suggests the existence of con- 
served regulatory mechanism for some AGOs in S. 
miltiorrhiza and Arabidopsis. 

It is generally considered that miRNA expression is 
negatively correlated with that of targeted mRNAs [50]. 
However, there are many exceptions. For instance, 
miR160, miR164 and miR172 showed a positive correl- 
ation with their targets in rice [51]. Comparing the expres- 
sion of S. miltiorrhiza miRNAs and their targets showed 
that no simple linear correlations existed between miR168 
and SmAGOl and between miR403 and SmAG02 



(Figures 4 and 7). It could be a consequence of the feed- 
back regulation of miR168/AG01 and miR403/AGO2 
[21,22]. It is also possible that these miRNAs and targets 
are regulated by other unidentified factors associated with 
a more complex regulation [52,53] . 

Although we cloned cDNA fragments with the 5 '-end 
mapped to the miR169* complementary region of 
SmAGO 5, they were probably not the products of 
miR169*-directed cleavage. First, the penalty scores for 
miR169*:«SmAG05 were at least 9.5, suggesting low com- 
plementarity between miR169* and SmAGOS. Second, 
the positions mapped were not located between the 10 
and 11 nucleotides from the 5 '-end of miR169*, which 
was atypical for miRNA-directed cleavage [54]. 

S. miltiorrhiza miR167*, miR396* and miR530 were 
computationally predicted to target SmAG08, SmAG03 
and SmAGOl for cleavage, but no 5-RACE products 
were obtained for these SmAGOs. It indicates that these 
miRNAs are probably not involved in the regulation of 
SmAGOs. However, we can not rule out the possibility 
that some of the miRNAs regulate SmAGO mRNA levels 
in tissues unanalyzed or at specific developmental stages 
of S. miltiorrhiza or in cells undergoing certain environ- 
mental stresses. It is also possible that some of the 
miRNAs interact with SmAGO mRNAs for translational 
repression. Further analyzing the interaction between 
miRNAs and SmAGOs through genetic transformation 
may give us a clearer picture about the regulatory mech- 
anism of SmAGOs. 

Conclusions 

The decoding of S. miltiorrhiza genome allowed us to per- 
form a genome-wide prediction of SmAGO genes. The pre- 
diction was further confirmed by full-length cDNA cloning, 
which resulted in identification of the first set of full-length 
AGO cDNAs in a plant species. Through a comprehensive 
approach, combining sequence feature, gene structure and 
conserved domain comparison of SmAGOs and AtAGOs; 
phylogenetic analysis of AGO proteins in S. miltiorrhiza, 
Arabidopsis and rice; and expression of SmAGO genes S. 
miltiorrhiza, we characterized the newly identified 10 
SmAGOs. The results suggest the conservation and diver- 
sity of sequence features and gene functions among AGOs 
from different plant species. Analysis of posttranscriptional 
regulation of AGO genes suggested the existence of con- 
served AG01-miR168 and AGO2-miR403 regulatory mod- 
ules in S. miltiorrhiza and Arabidopsis. These results will 
help to open a window for understanding gene silencing 
networks in the well-known TCM plant, S. miltiorrhiza. 

Methods 

Plant materials 

S. miltiorrhiza Bunge (line 993) with whole genome se- 
quences available was grown in a field nursery. Mature 
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flower buds blooming soon, mature and healthy leaves, 
young stems and roots in about 0.5 cm diameter were 
collected from 2-year-old plants on Aug. 15, 2012. Sam- 
ples from three plants were pooled and stored in liquid 
nitrogen until use. 

Prediction of SmAGO genes 

Arabidopsis and rice AGO protein sequences were 
downloaded from GenBank (http://www.ncbi.nlm.nih. 
gov/protein) and used to search for homologues against 
the current assembly of the S. miltiorrhiza genome 
(Chen et al., unpublished) using the tBLASTn algorithm 
[27]. An e-value cut-off of 10" 10 was applied to the 
homologue recognition. All retrieved sequences were 
used for gene prediction on the Genscan web server 
(http://genes.mit.edu/GENSCAN.html) [28]. The pre- 
dicted gene models were further examined and corrected 
manually by comparison with AGO genes identified 
from other plant species using the BLASTx algorithm 
(http://www.ncbi.nlm.nih.gov/BLAST) [27] . 

Cloning of SmAGO genes 

Total RNA was extracted from the root of S. miltiorrhiza 
using Trizol reagent (invitrogen). mRNA was purified 
using the oligotex mRNA mini kit (invitrogen). RNA 
ligase-mediated rapid amplification of 5' cDNA ends (5'- 
RACE) and 3' cDNA ends (3 '-RACE) was performed on 
mRNA using the GeneRacer kit (Invitrogen). PCR amplifi- 
cations were performed using the GeneRacer primer and 
the nesting gene-specific primers (see Additional files 3 
and 4) under the following conditions: predenaturation at 
94°C for 2 min, 5 cycles of amplification at 94°C for 30 s 
and 72°C for 1 min, 5 cycles of amplification at 94°C for 
30 s and 70°C for 1 min, 25 cycles of amplification at 94°C 
for 30 s, 56°C for 30 s and 72°C for 2 min, followed by a 
final extension at 72°C for 15 min. Nested PCR amplifica- 
tions were carried out using the GeneRacer nested primer 
and the nested gene-specific primers (see Additional files 
3 and 4) under the following conditions: predenaturation 
at 94°C for 2 min, 30 cycles of amplification at 94°C for 
30 s, 58°C for 30 s and 72°C for 2 min, followed by a final 
extension at 72°C for 15 min. PCR products were purified, 
cloned and sequenced. 

Based on the obtained 5 ' and 3 ' cDNA sequence, gene- 
specific forward and reverse primers (see Additional file 5) 
were designed for amplification of full-length SmAGO 
cDNA. PCR amplifications were performed under the 
following conditions: predenaturation at 94°C for 2 min, 
30 cycles of amplification at 94°C for 30 s, 56°C for 30 s 
and 72°C for 3 min, followed by a final extension at 72°C 
for 15 min. PCR products were gel-purified, cloned 
and sequenced. 



Bioinformatic analysis and phylogenetic tree construction 

The molecular weight (MW) and theoretical isoelectric 
point (p7) were predicted using the Compute pI/MW 
tool on the ExPASy server (http://web.expasy.org/ 
compute_pi/). Intron/exon structures were analyzed 
manually based on genomic DNA sequences and the 
cloned cDNA sequences. Conserved domains were ana- 
lyzed by search the deduced amino acid sequence of 
SmAGOs against the NCBI conserved domain (http:// 
www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi). The con- 
served residues were analyzed by alignment of amino acid 
sequences using DNAMAN and then checked manually. 
For phylogenetic tree construction, amino acid sequences 
of AGOs from S. miltiorrhiza, Arabidopsis and rice 
were first aligned using ClustalW version 1.83 and 
then constructed using MEGA version 4.0 by the 
neighbor-joining method with bootstrap to be 1000 
replicates [55,56]. 

Quantitative real-time reverse transcription-PCR (qRT-PCR) 

Total RNA was extracted from plant tissues using the 
plant total RNA extraction kit (BioTeke) and then 
treated with RNase-free DNase (Promega) to remove 
genomic DNA contamination. RNA integrity was ana- 
lyzed on a 1% agarose gel. RNA quantity was determined 
using a NanoDrop 2000C Spectrophotometer (Thermo 
Scientific). Reverse transcription was performed on 1 ug 
total RNA by 200 U Superscript III reverse transcriptase 
(Invitrogen) in a 20 ul volume. The resulting cDNA was 
diluted to 200 ul with sterile water. qPCR was carried 
out in triplicates using the BIO-RAD CFX system (Bio- 
Rad). Gene-specific primers were listed in Additional 
file 6. SmUBQIO was used as a reference as previously 
described [26]. PCR was carried out in a 20 ul volume 
containing 2 ul diluted cDNA, 250 nM forward primer, 
250 nM reverse primer, and 1 x SYBR Premix Ex Taq II 
(TaKaRa) using the following conditions: predenatura- 
tion at 95°C for 30 s, 40 cycles of amplification at 95°C 
for 5 s, 60°C for 18 s and 72°C for 15 s. The results from 
gene-specific amplification were analyzed using the com- 
parative Cq method which uses an arithmetic formula, 
2-AACq, to achieve results for relative quantification 
[57]. Cq represents the threshold cycle. 

The levels of miR168a/b and miR403 were analyzed 
using the miRNA-specific poly(T) adaptor RT-PCR 
method [40]. Briefly, 1 ug DNasel-treated total RNA 
was polyadenylated using the Poly(A) Tailing kit 
(Ambion) as described previously [40] and then 
reverse-transcribed into single-strand cDNA. qRT- 
PCRs were performed in triplicates as described [58]. 
S. miltiorrhiza 5.8S rRNA was used as an endogenous 
reference. Primers used for miRNA quantification 
were listed in Additional file 7. 
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Identification of S. miltiorrhiza miRNAs with perfect or 
near-perfect complementarity to SmAGOs 

Known plant miRNA sequences were downloaded from 
miRBase (release 19, http://www.mirbase.org/) [35]. Target 
search of known plant miRNAs was performed against ten 
full-length SmAGO cDNA sequences using psRNATarget 
[36]. The maximum expectations of 3.5 and the target 
accessibility-allowed maximum energy to unpair the target 
site of 50 were applied. The identified plant miRNA se- 
quences were then aligned with the current assembly of 
the S. miltiorrhiza genome (Chen et al, unpublished) 
using the SOAP2 program with no more than 2 mis- 
matches allowed [37]. Hairpin structures were predicted 
using the mfold program [38]. Criteria described by [59] 
were applied to annotate S. miltiorrhiza miRNAs. 

Mapping of SmAGO cleavage sites 

SmAGO cleavage sites were mapped using the modified 
RNA ligase-mediated rapid amplification of 5' cDNAs 
method as described [39]. PCRs were carried out on 
mRNA isolated from S. miltiorrhiza roots and pooled 
samples containing flowers, leaves, stem and roots using 
the GeneRacer 5' primer and the nesting gene-specific 
primers (see Additional file 8). Nested PCRs were 
performed using the GeneRacer 5' nested primer and 
the nested gene-specific primers (see Additional file 8). 

Additional files 
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